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Abstract 

The supergravity dual of N regular and M fractional Dl-branes on the cone over the 
Einstein manifold Q 1,1,1 has a naked singularity in the infrared. The supergravity dual 
of N regular and M fractional D3-branes on the conifold also has such a singularity. 
Buchel suggested and Gubser et al. have shown that in the D3-brane case, the naked 
singularity is cloaked by a horizon at a sufficiently high temperature. In this paper we 
derive the system of second-order differential equations necessary to find such a solution 
for Q 1,1,1 . We also find solutions to this system in perturbation theory that is valid when 
the Hawking temperature of the horizon is very high. 
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1 Introduction 



The AdS/CFT correspondence JT|, |], Q has produced a wealth of new information about 
strongly coupled conformal gauge theories. Considerable effort has also been invested 
into extending it to non- conformal theories. One recent development is the fascinating 
story that has emerged surrounding a certain four dimensional M = 1 supersymmetric 
SU (N) x SU (N + M) gauge theory 0, |5| . The goal of the present work is to attempt to 
retell at least part of the same story || for a M = 2 supersymmetric SU(N) x SU(N) x 
SU(N + M) gauge theory living in two dimensions. 

We begin by summarizing the story surrounding the M = 1 SU(N) x SU (N + M) 
gauge theory. The theory may be realized by adding M fractional D3-branes (wrapped 
D5-branes) to N regular D3-branes at the apex of the conifold, which is defined by 
the constraint Y^=i z f = in C 4 0. For M = 0, this gauge theory reduces to the 
superconformal theory dual to the AdS§ x T 1,1 background of type IIB string theory 
01- 

In the supergravity dual, the M fractional branes correspond to M units of RR 3-form 
flux through the 3-cycle of the compact space T 1,1 . This flux changes the background 
and introduces the logarithmic running of J S 2 B2, which is related to the running of 
field theoretic couplings ||. In turn, this running causes the RR 5- form flux, which 
corresponds to the number of ordinary D3-branes, to grow logarithmically with the radius 
0, due to the equation dF 5 = H 3 A F 3 . In ||, this behavior was attributed to a cascade 
of Seiberg dualities in the dual gauge theory. 

While the Klebanov-Tseytlin (KT) solution || is smooth in the UV (for large p), it 
has a naked singularity in the IR. Two complementary ways have been found of removing 
the singularity, and it is with the removal of the singularity that the story becomes very 
interesting. In J4j, Klebanov and Strassler (KS) proposed to deform the conifold, i.e. 
to replace the constraint with J2i=i z \ — ^ ■ The resulting solution, a warped deformed 
conifold, is perfectly non-singular and without a horizon in the IR, while it asymptotically 
approaches the KT solution in the UV. The mechanism that removes the naked 
singularity is related to the breaking of the chiral symmetry in the dual SU(N) x SU (N + 
M) gauge theory. The Z 2 m chiral symmetry, which may be approximated by U(l) for 
large M, is realized geometrically e . The deformation of the conifold breaks 

it down to Z 2 : z$ — > — z< although supersymmetry is preserved. 

The second mechanism for removing the singularity from the KT solution was pro- 
posed by Buchel [10| and later worked out in detail [|ll|, ||. It was suggested that a 
non-extremal and hence supersymmetry breaking generalization of the KT solution, with 
unbroken U(l) symmetry, may have a regular Schwarzschild horizon "cloaking" the naked 
singularity. The dual field theory interpretation is restoration of chiral symmetry above 



some critical temperature T c . In [11, H], the authors were able to show that at least at 



high temperatures, where the differential equations could be analyzed through perturba- 
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tion theory in the number of fractional D3-branes, a well behaved supergravity solution 
exists with restored U(l) symmetry which in the IR involves a regular Schwarzschild 
horizon but which in the UV approaches the asymptotic KT geometry. 

These two methods of removing the singularity from the KT solution form an at- 
tractive and consistent picture for the gauge theory dual. At high temperature, above 
T c , the chiral U(l) symmetry is present. As we lower the temperature, the horizon dis- 
tance shrinks until we reach the critical temperature T c where the horizon can no longer 
"shield" the singularity. At this point, a phase transition occurs and a KS type solution 
0| becomes preferred. The chiral symmetry is broken. 

In Jl2|], generalizations of the KT solution were found involving fractional M2-branes 
and fractional Dp-branes, p = 0, 1,2,4. Like in the KT solution, the transverse space is 
conical and moreover these generalizations are typically smooth in the UV but possess 
a naked singularity as the radius of the cone shrinks below some critical value. This 
singularity renders the dual gauge theory poorly defined in the IR. A detailed under- 
standing of the gauge/gravity correspondence for these fractional branes could shed light 
on gauge theories in dimensions other than four, and so removing the IR singularities 
of these generalizations is an important challenge. In the case of fractional D2-branes, 
a KS type @] solution has already been found |R| that resolves the singularity while 
preserving supersymmetry. The IR limit p — > of this fractional D2-brane solution is 
thought to correspond to confinement in the dual three dimensional gauge theory. 

It is natural to wonder if the mechanism for singularity resolution at high Hawking 
temperature in the fractional D3-brane system will work for fractional Dp-branes with p 
other than 3.Q We shall focus on the case p — 1, the fractional D-string solution. This 
solution consists of a warped product of R 1 ' 1 flat space-time directions and a Ricci flat, 
eight dimensional cone Y 8 . The base of the cone is then a seven dimensional Einstein 
space, which in this paper will be Q 1,1,1 [14j]. The manifold Q 1 ' 1 ' 1 can be described as a 
coset: 

Xil)1 = SU(2f 
V U(l) 2 ' 

The ordinary D-strings fill the space-time dimensions. D3-branes are wrapped over one 
of the 2-cycles of Q 1 ' 1 ' , and the remaining two dimensions of these D3-branes fill R 1,1 , 
creating the fractional D-strings. In the dual supergravity description, the fractional 
D-strings correspond to turning on an RR five form flux, F 5 , which pierces one of the 
five-cycles of Q 1,1,1 . Meanwhile, the ordinary D-strings correspond to an electric form 
flux F 3 in the R 1 ' 1 and radial directions. 

Our non-extremal generalization of the fractional Dl-brane solution is well-behaved 
and singularity free at high Hawking temperature. Moreover, our ansatz preserves the 



1 A related finite temperature solution was discussed in |15| , which used the wrapped brane solution 
of |16|] to identify appropriate boundary conditions in the IR. 
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U(l) symmetry of Q 1,1,x . Hopefully, a U(l) breaking deformation of the cone over Q 1,1,1 , 
which eliminates the singularity of the KT type fractional D-string solution, will someday 
be found. It seems reasonable to conjecture that there is some chiral symmetry breaking 
phase transition, just as in the fractional D3-brane case f5|. 

Below the critical temperature T c , this hypothetical deformed cone over Q 1 ' 1,1 would 
be preferred. Above T c , the ansatz we present below describes a fractional D-string 
solution where the naked singularity is cloaked and the U(l) symmetry restored. 

We begin by presenting our ansatz for this U(l) symmetry preserving fractional D- 
string. Next, we develop a perturbation theory that is valid at high temperature, and 
thereby show that at sufficiently high temperature, the naked singularity of the corre- 
sponding KT type fractional D-string is "shielded" by an event horizon. 

2 Non- Extremal Generalization of the Fractional Dl- 
Brane Ansatz 

We start with an ansatz for the non-extremal fractional Dl-branes. Our strategy is 
similar to that used in finding the non-extremal generalization of the fractional D3- 
branes Jllfl . We will add additional warping functions to the metric that preserve the 
underlying U(l) x SU(2) 3 symmetry of Q 1,1,1 . At the same time, we will leave unchanged 
the geometric dependence of the RR and NSNS field strengths on the two- and five-cycles 
of Q 1 ' 1 ' 1 . 

The general ansatz for a 10-d Einstein-frame metric consistent with the underlying 
symmetries of Q 1 ' 1 ' 1 involves five functions x, y, z, w\, and W2 of a radial coordinate u 



where 



and 



dsi 0E = e iz (e~ 2x dX 2 + e 2x dXf) + e~ z dsl (1) 

ds 2 % = e Uy du 2 + e 2y (dM 7 ) 2 , (2) 
{dM 7 f = e" 12 -4 + e 2 -" 2 - J2(gl + gfo + e 2 ^ 4 ^ + g^) (3) 

i=l 

1 3 1 1 



Here X is the Euclidean time and Xi is the longitudinal D-string direction. 



2 Although phase transitions are usually forbidden for 1+1 dimensional systems, they can occur for 
gauge theories in the limit of infinite N jL7). This phase transition is a large N phenomenon; for 
finite N we would expect a smooth crossover between the symmetric high-temperature phase and the 
low-temperature phase. 
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This metric can be brought into a more familiar D-string form 

dp 2 



ds 2 10E = h{p)-*l A [A{p)dXl + dXl] + h{pf' A 



B(P) 



+ p 2 (dM 7 



(5) 



with the redefinitions 

h = e -4,-8x/3 ^ p = e y +X /3 ? A = ? ^+2^/3^2 = ^-1^2 (g) 

When wi = u? 2 = and e 6y = p 6 = ^, the transverse 8-d space is the cone over 
Q 1 ' 1,1 = M 7 . Small u thus corresponds to large distances (where we shall assume that h, 
A, B — > 1, and p — > oo) and vice versa. 

The function w\ squashes the U(l) fiber of Q 1 ' 1 ' 1 relative to the three two spheres 
while the function w 2 squashes one of the two spheres relative to the other two. For 
comparison, the non-extremal generalization of the KT solution, which involved the 5-d 
Einstein manifold T 1,1 , made use of only one squashing function w, not two. This w, 
which squashed the U(l) fiber of T 1,1 relative to rest of the manifold, roughly corresponds 
to our w\. The most general volume preserving squashing of Q 1 ' ' consistent with the 
symmetries would involve two more functions Wi and 103. Note, however, from (|8|), that 
the harmonic two- forms on Q 1 ' 1 ' 1 involve linear combinations of the SU (2)'s. In order 
to keep one of these harmonic two- forms u independent, we may only add one more 
squashing function w 2 . The additional Einstein equations in the Q 1 ' 1 ' 1 directions give 
additional constraints when compared to the five dimensional T 1,1 case; fortunately w± 
and W2 give us the freedom necessary to satisfy the constraints. Ricci-flat 8-d spaces with 
nontrivial Wi correspond to the resolutions of Q 1 ' 1 ' 1 considered by |L8 . 



The extremal D-string solution and the more general fractional D-string solution on 
Yg have x = Wi = 0. Adding a non-constant x(u) drives the non-extremality. For 
example, the non-extremal version of a D-string on the cone over Q 1,1 ' 1 (w = 0) has 
x = au, e~ 4x = A = 1 - |%, e" 42 " 8 ^ 3 = h= l + ^, p = e y+x ^. Our aim will be 
to understand how switching on the non-extremality (x = au) changes the extremal 
fractional D-string solution. 

Our ansatz for the p-form fields is dictated by the geometry and thus is exactly the 
same as in the extremal fractional D-string case ||12|| : 



K{u)e^d 2 x A du , (7) 



B 2 = f(u)u 2 , u 2 = (g 6l A - g 02 A g <t>2 )/\ / 2 , (8) 

F 5 = P(uj 5 + *u 5 ) , u 5 = -sy, Au 2 A g 03 A g^ . (9) 

Moreover, the dilaton $ is assumed to be a function of the radial variable u only. The M 
fractional D-strings (wrapped D3-branes) thus correspond M units of flux through the 
five-cycle of Q 1 ' 1 ' 1 , and P ~ g s M. The function K{u) corresponds roughly to the flux of 
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ordinary D-strings through the compact space Q 1,1 ' 1 itself. The equation of motion for 
F 3 , d*e*F 3 = F 5 A H 3 , implies 

K(u) = Q + Pf(u) . (10) 

In what follows, we shall use this ansatz to reduce the type IIB supergravity equations 
of motion to a system of nonlinear, coupled ordinary differential equations describing the 
radial evolution of x, y, z, Wi, w 2 , K, and $. 



3 Derivation of the Equations of Motion 

We have seven warping functions and hence will require a system of seven ordinary differ- 
ential equations. From analogy with the non-extremal generalization of the KT solution 



11], we also expect a zero energy constraint, giving eight equations total. Consideration 



of the p-form field strengths yields two nontrivial equations of motion, one for the dilaton 
and one for H 3 = dB 2 . The H 3 equation of motion, d*e~®H 3 = — F 5 A F 3 , reduces to the 
ordinary differential equation 

while the dilaton equation of motion, 2 d*d<& = —e~®H 3 A *H 3 + e®F 3 A *F 3 , reduces to 

K' 2 

2Q" — _ e -$+2z-4y-4w 1 +4w2 ^2 g 6.z-$ /-j^N 

p2 

Einstein's equations are Rmn = Tmn where Rmn is the Ricci curvature and 

Tmn = -9m^9 n ^ + —F MP q RS F n pqrs 
2 9o 

+^( e_<1> HmpqHn PQ + F M pqF n pq ) 

-^G MN (e-*H PQR H p Q R + e*F PQR F p Q R ) . (13) 

In order to write down these equations in a convenient form, we will work in an orthonor- 
mal frame basis: 

e° = e^ z - x dX , e 1 = e^ +x dX t , e u = e~^ z+7y du , (14) 



e * = e 



-\Z +V+Wl - W2g ^ ) ^ = e -l Z+ y +Wl - W2g ^ ^ i= l i2 . (16) 



In this basis, Einstein's equations are diagonal. Moreover, the equations corresponding to 
Re^, Re 2 e 2 , R(j>i<t>ii an d R<j> 2 ij> 2 are identical and similarly for the equations corresponding 
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to Re 3 e 3 and Rfafo. Thus, we are left with six relations. Our strategy will be to put 
aside R uu at first and use the remaining five relations along with the two field strength 
equations to derive a second order, nonlinear system in the seven warping functions. At 
the end, we will find that the R uu relation provides a zero energy constraint, analogous 



to the one found in [ |TT|j . 

From the field strengths, it is relatively easy to see that T 00 = T\\. However, R 00 = 
— e z ~ My (^z"+x"), while R\\ = —e z ~ 1Ay {^z" — x"). Hence, the first two Einstein equations 
allow us to solve for x(u) exactly: 

x" — , x = au , a > . (17) 



In the case of non-extremal fractional D3-branes [HJ, the same behavior was found for 
this x function, and the factor a was identified with the degree of non-extremality. 

Having solved for x, we can use either of the first two Einstein equations to find an 
equation for z: 

I2z" = 2P 2 e 4z+4:y+Awi ~ Aw2 + 3K 2 e 6z ~® + — - c 2z ~ 4 ^~ 4w i+ 4w 2~^ (ig) 



Comparing (|T2]) and ( pT8[) with the corresponding equations (3.18) and (3.19) in 
there is some interesting similarity. Both here and in [|ll[], the equations involve sums of 
the same three terms proportional to P 2 , K 2 , and (K'/P) 2 . In fact, there exists a linear 
transformation of (|Tj|) and (|18D , 

1 3 _1 

$ = U> n - 2z n , $ n = U> + z , (19) 



such that the new differential equations (|23|) and (p4[) correspond almost precisely with 



those in JTTj] . To take advantage of the calculations in |Tl| and ||, we shall use the 
transformed variables z n and $ n in what follows. Note that if = 0, then z" = — <J>"/2. 
Moreover, in the extremal case x = 0, z = — $/2 and h~ 3 ^ = e 3z = h^ x l 2 e^l 2 . Hence, 
we see that the Einstein frame metric ([|) in this particular case can be obtained from 
the string frame metric through the usual procedure of multiplying by a factor of e^*/ 2 . 

To derive similar equations for y, Wi, and W2, we need to take linear combinations of 
the Einstein equations involving 

Rw = e z - Uy (^z" - y" + 6w'() + 2e z ~ 2y ■( e -^i-8« 2 + 2e -ie« 1 +4 tt2 ) ? 
R 6i9i = e z - Uy {^z" - y" - w'[ + w'l) + 2e z - 2y (-e- Ww ' +4w2 + Ae~ 2w ' +2w2 ) , 
R e303 = e z - Uy {-z" - y" - w'l - 2wl) + 2e z ~ 2y (_ e -i^i-8^ + 4e ~2 Wl -4 W ^ _ 
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We leave it to the reader to derive the field strength contributions ([13]) to Einstein's 
equations and give only the results 



7y" = + 2e 12y (- e - Wwi - Sw2 - 2e~ mwi+4w2 + 8 e - 2wi - 4w2 + 16 e - 2wi+2w2 ) , (20) 

Q 

Yw" = $" + —e 12y (_ e -16u, i- 8,u 2 _ 2g- 16tu i+ 4tu 2 _|_ e -2tui-4tu 2 _|_ 2 e ~ 2w i+ 2w 2^ ^21) 

2^2 = - l e 12^( e -16^i-8^ 2 _ e -16^ 1+ 4«, 2 _ 4e -2«, 1 -4«, 2 + 4e -2«, 1+ 2™ 2 ) ^ (22) 

(K>) 2 

e 42 n -4j/-$„-4«)i+4u)2 ^23) 

g^" _ 2K 2 e 8Zn + p 2 e 42 "+ 4 ?'+*™+ 4w i- 4,i '2 _|_ ) c 4:Zn-4y-&„-4wi+4w2 ^4) 

^ e iz n -4y-'i> n ~Aw 1 +Aw2 j£f\l _ p 2 J{ e Sz ™ (25) 

As yet, we have ignored the R uu Einstein equation. For our metric 

Run = e- 1 ^ (\z" - Q(z') 2 - 7y" + A2(y') 2 - 2{x'f - A2(w[) 2 - 12(vQ 2 ) (26) 



which, using (p0|) , (^4|), (p3|) , (|19D and (p~7|) , yields the zero energy constraint 
8%') 2 - 16«) 2 - 3($^) 2 - 84K) 2 - 24K) 2 + 6< + K 2 e 8 -- 

4g 12 f(_ e _16w 'i _8w '2 _ 2g- 16 ' !u i+ 4u '2 _|_ g e -2iDi-4«;2 _|_ j^g e -2toi+2ii;2^ _ 4fl 2 (27) 

In later sections, it will be important to keep track of the dimensions of the various 
parameters involved. From the form of the metric (H) it is natural to require that e v and 
m -1 / 6 should have dimension of length, while x, z, w be dimensionless. Since we have set 
the 10-d gravitational constant to be 1 (i.e. we measure the scales in terms of the 10-d 
"Planck scale" Lp ~ (g^c/ 2 ) 1 / 4 ), then from fl2"T|), we conclude that K and Q in flip]) have 
dimension (length) 6 while P has dimension (length) 4 and / has dimension (length) 2 . It 
is easy to restore the dependence on the Planck length by rescaling Q — > LpQ, P — > 
LpP, etc. To restore the dependence on the string coupling one should further rescale 
P 2 — > g s -P 2 . At the end, Q ~ g s a' 3 N, P ~ g s a' 2 M } where iV and M are the number of 
ordinary and fractional D-strings, respectively. 



4 Three Simple Solutions 

In addition to the extremal D-string solution, there are three other relatively simple 
solutions to the system of equations (p0|)-([25|), ([271). Some of of these solutions were 
discussed briefly in the Introduction and in Section 2. There is the analog of the KT 
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solution for fractional D-strings, what we will call the extremal fractional D-string so- 
lution. Next, there is the non-extremal ordinary D-string solution. These two solutions 
will be very important when we try to find a non-extremal fractional D-string solution 
through perturbing in the number of P of fractional D-strings. We will find a solution 
which interpolates between the extremal fractional D-string solution in the UV and the 
non-extremal ordinary D-string solution in the infrared. 

The third solution is more of a mathematical curiosity. It is the analog of the sin- 
gular, non-extremal D3-brane solution found in 1 10(1 - This third solution, although non- 



extremal, is singular because it has a naked singularity in the far infrared. Although it 
approaches the extremal fractional D-string solution as the non-extremality parameter 
a — > 0, it does not approach the non-extremal ordinary D-string solution as the number 
of fractional D-strings P — > 0. 

4.1 Singular Non-Extremal Fractional D-Strings 

As a first attempt at finding a non-extremal solution, we might be tempted to try to 
preserve the geometry of the Q > ' base space by turning off the squashing functions wf, 
this approach will lead us to the D-string analog of the Buchel solution. Our motivation 
is two-fold. First, because this solution is singular, we will see the necessity of squashing 
the Q ' ' . Second, to obtain the extremal fractional D-strings in the next subsection, we 
need only take the limit in which the non-extremality parameter a — > 0. 

So let us suppose that the Q 1 ' 1 ' 1 is not squashed, i.e. Wi = w 2 = 0. Then the 
equations (|20|)-(p5|), (|27|) simplify dramatically. From (|22|) , we have Q' n = p, where p is 
a constant, so fl23|) becomes 

f = _p e %+*n . (28) 
From (|25|), it is straightforward to see that z n must then also satisfy a first order equation: 

{e-^)' = Q + Pf . (29) 

Then ( pU|) reduces to y" = 6e 12y . The zero-energy constraint (^7]) sets one of the integra- 



tion constants of this differential equation: 



y ' = -V6 2 + e 12 y , 84b 2 = 4a 2 + 3p 2 , (30) 

which integrates to give e Sy = sin ^ 6bu ■ Without loss of generality, we may choose b > 0. 
The differential equation for / ( p8|) becomes 

/' = -p^- a) " (t3^)' 3 ■ (3i) 

Once / is known, z n is easily found by integrating fl29"|). 
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The precise analog with the Buchel solution is found by taking p = — 2a/ 3 which in 
turn implies that 6 = 2a/3\/7 ([301) , but we shall keep p and 6 general in the discussion 
that follows. To demonstrate the pathological nature of these solutions, let us take a 
look at the Ricci scalar: 

R = I e *.-i4irtW4 ^p2 e ^+4z n+P u - I ((Q + P/) 2 e 82 " - Qp(Q + Pf)e 4z " - 9p 2 )"j (32) 

In analogy with the behavior of the Buchel solution for D3-branes, we expect this solution 
will exhibit divergences at large u, so we develop asymptotics for / and z n far in the IR: 

/ = /„ - P ^ 2 ^ z {v ~ Ab)u + 0(e ip - 16b)u ) (33) 

and hence 

p2(2h\2/3 

e~ 4z - =Ci + u(Q + Pfo) - ? Auv> e<Jl> ~ U + 0(e^- 16b >) (34) 

[p — 46) z 

where Ci is some integration constant. The asymptotics are clearly very sensitive to the 
relative magnitudes of p and b. Hence we will consider the two cases, p — 46 > and 
p — 46 < separately. 

In the case p — 46 > 0, the exponential term in fl34|) is dominant and causes e _4z ™ 
to be large and negative at large u. The subsequent exponential terms contribute very 
little because of the constraint fl3~0|): in particular, in order for a to be real, \p\ < \/286 
and hence p — 166 < 0. The function e~ 4zn is continuous in the region < u < oo, and 
note that near u — 0, e~ 42:n ~ C 2 + it(Q + P/o) + 0(u 4//3 ). Hence, close to u — 0, e~ AZn is 
nonnegative (in order to have well behaved asymptotics at large distances) and growing. 
From these facts, it is clear that e~ 4Zn becomes zero at some finite u cr and the Ricci 
scalar blows up. Defining the horizon as the locus where G 00 = exp(3^„ — 2au + 3pw/4) 
vanishes, a horizon will only occur, if at all, in the limit u — ► oo. Hence, the singularity 
is naked and the solution is pathological. 

For the case p — 46 < 0, the constant and linear terms in ( |3"3"D and (K3j) dominate at 
large u. Hence, / ~ fo and e ™ ~ u. Also, e~ v ~ e bu . As was noted previously, in order 
for a to be real, \/28b > \p\. Hence, the e v terms will dominate the Ricci scalar and send 
it off to infinity at large u, even if we consider the limit P — ► 0. In other words, this 
solution has a naked singularity in the IR, even in the absence of fractional D-strings. 
The best we can do is to set the relative magnitudes of the a and p so that there is also a 
horizon at u — oo. These considerations show that leaving the Q 1 ' 1,1 unsquashed leads to 
a naked singularity, so to find a well behaved non-extremal fractional D-string solution, 
we will have to look elsewhere. 
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4.2 The Extremal Fractional D-string Solution 



As was mentioned above, the limit a, p — > of the singular non-extremal D-string solution 
recovers the extremal fractional D-string solution. Let us take this limit: 

*»-0, e 6 ^-^, f^f -^ u y\ (35) 
bu 2 

3 • fi 1 / 3 P 2 

^ Ci + {Q + Pfo)u M 4/ 3 (36) 

o 

This solution is well behaved in the UV, in the limit u —>■ 0. If C\ — 1, then the metric 
(§) becomes asymptotically flat at large distance u — > 0. However, in what follows, it will 
be more convenient to choose C\ = 0, thus eliminating the asymptotically flat regime, 
"zooming in" on the low-energy dynamics of the dual gauge theory. As was noted in |T^| . 
the solution gives rise to a naked sing ularity in the infrared at u\f = (Q+P/ )4-6 2/3 /9P 2 . 
In analogy with , we shall define Lp as the value u„V 6 at which this solution develops 
a singularity. 



4.3 The Non-Extremal Ordinary D-String Solution 

We are searching for a solution without fractional D-strings, so we set / = P = 0. 
Moreover, we know from the literature |19[ that this solution does not require the extra 
degrees of freedom provided by w± and wi so we set w\ = w<i = 0. The system of equations 
(p0|) — (|25|) , (|2T| ) simplifies substantially in this case and becomes easily integrable: 

y" = Qe 12 y, < = ^QV-, x" = , < = . (37) 
Integrating these equations once, we find 

x' = a, K=P, y' 2 = b* + e^, < 2 = c 2 + ^e 8 -. (38) 

lb 

From the constraint fl27D , we find that the integration constants must obey a 2 — 21b 2 + 
4c 2 + \p 2 = 0. Integrating the equations for y and z n , we find that 

e 6y = b , e AZn = 4c . (39) 

sinh 66m ' Q sinh 4c(u + k) 



Cast in a more familiar form, 



P 6 = * 6y+2x = , „_ 12bu , (40) 



2^g-2(36-a)u 

1-e- 

A(u) = e~ 4x = e~ 4au , (41) 
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h( U ) = e -^-*n-8x/3 = Q_ e - pu -8au/3 ^ ^ + fc ) _ ^ 

As has become customary in gauge/gravity duality we set the boundary condition that 
h approach zero as u — > 0. In other words k = 0. This boundary condition removes the 
asymptotically flat region so that we "zoom in" on the low-energy dynamics of the dual 
gauge theory. 

To match the standard non-extremal D-string solution, we require that 36 = 2c = a, 
and p = — 2a/3. This choice satisfies the zero-energy constraint and (f40|)-(|42"|) become 

, A{p) = B{p) = l-— R 



KP) = £e , ^(p) = = 1 - — 6 . (43) 



In other words, we recover the standard non-extremal D-string metric. To summarize, 
the non-extremal ordinary D-string solution consists of Wi = and 

$ n = -2au/3 , e 4x = e 4au , e~ 6y = 3a" 1 sinh2aw , e" 42 " = £ sinh2au . (44) 

5 Asymptotics of the Regular Non-Extremal Frac- 
tional D-String 

We have not succeeded in finding an analytic solution to the system of differential equa- 
tions (p0|)-(|25|), (f27|). To proceed further, we can either integrate the equations numer- 
ically or seek an approximate solution in perturbation theory. In any case we must be 
certain that our solution satisfies the correct boundary conditions in the short- distance 
(u — * oo) and long-distance (u — ■> 0) limits, where we understand the physics. The 
procedure we will follow is very similar to that in ||. 

For P —>■ we must obtain the black D-string solution, which has a regular Schwarzschild 
horizon. If the horizon is preserved as we add fractional Dl-branes, we expect the fol- 



lowing asymptotics to hold as u — > oo ((44|): 

x = au , y — > —aw/3 + y* , z n — > —au/2 + , 

Wi -> Wi* , $ n -> -2aw/3 + , K -*K m . (45) 
The metric for w — > oo in the w — Xo directions is given by 

ds 2 = e- 4au+3z * + ^*dX 2 + e- 4au - z *-^* +Uy *du 2 (46) 

so that with the natural near- horizon variable U = e~ 2au , the usual procedure of choosing 
periodicity of the Euclidean time X to avoid a conical singularity fixes the Hawking 
temperature T#: 
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Tjj = l e ^-7y*+**/2 . (47) 

7T 

At large distances (u — > 0) the non-extremal solution should approach the extremal 
solution (p5|), (p6|), i.e. we require that 

it — > : x, »0, y — »■ — log 6m . (48) 

6 



(Note that this behavior is also in agreement with the small u asymptotics (f44|) of the 
regular non-extremal D-string solution.) The behaviors of the effective D-string charge 
and of the warp factor at small u are 

«l/3p2 / o \ q . Rl/3p2 

K{u) - °-f- (^L P 2 - u^) , - 6 -^-u(Lf - u^) . (49) 

We expect that the physics of this fractional D-string system should be very similar 
to that of the fractional D3-brane solution considered in [Bj. On the supergravity side, 
when T < T c the solution which preserves the U(l) symmetry of Q 1 ' 1 ' 1 is singular and 
one needs an appropriate deformation, perhaps of the KS-type [§|], that breaks this U(l) 
symmetry and removes the singularity. For T > T c we should be able to construct a 
nonsingular solution which preserves the U(l) chiral symmetry, and we will do this in 
perturbation theory in the next section. 

The flux corresponding to the number of D-strings is given by *F$ = K(u)e~ U2 Aws, 
so on the gauge theory side, we want to think of K(u)e~® , the effective D-string charge, 
as an effective number of unconfined color degrees of freedom. As we run the scale of 
the theory into the infrared (u — ► oo), the number of colors should decrease. Above 
the critical temperature, this number will be positive everywhere, but below the critical 
temperature, the effective number of colors will vanish at finite u. This behavior is the 
same as for fractional D3-branes. One potentially bothersome difference from the D3- 
brane case is the dependence of the D-string flux on the dilaton. However, things are all 
right because 



<K = e -f*n+2*„ (50) 



'sinh2a(ii + k) 

is decreasing for all positive u. Thus if K{u) decreases with increasing u, the flux will still 
decrease as well, provided that the fractional D-strings are a small enough perturbation 
that the variation of the dilaton is dominated by the presence of the ordinary D-strings. 
Moreover, K{u)e~® is well-defined for all values of u, in particular as u oo. Notice that 
in we have temporarily restored the integration constant k from (|3"9"|). In Section 4.3, 



we set k = to zoom in on the IR physics. This approximation introduces a singularity at 
u = (e - * diverges) but this is strictly an artifact of having removed the asymptotically 
flat region. 
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While the number of colors decreases as we move toward the horizon, from ( p0|) 
one can see that the string coupling e* increases. From [20], we expect that the string 
coupling can be expressed in terms of gauge theory parameters as e* ~ gy M N 1 / 2 / A 3 
where A sets the energy scale of the gauge theory, and N ~ Ke~® gives the number of 
ordinary D-strings at a given scale. Hence, e* ~ gyj^K 1 ^ 3 / A 2 . In the case where K is 
constant, i.e. there are no fractional D-strings, we expect the string coupling to become 
larger in the IR. Indeed, (|50"D holds exactly in this case. However, once we add fractional 
D-strings, K should decrease as u — > oo. Thus, in perturbation theory, we expect to see 
corrections to fl50[) that tend to decrease e* as u — > oo. 



6 Perturbation Theory in P 

One useful approach to constructing the required non-extremal fractional brane solution is 
to start with the non-extremal ordinary D-string solution (fP|), which is valid for P = 0, 
and find its deformation order by order in P 2 . A remarkable feature of perturbation 
theory in P 2 near the extremal (a = 0) D-string background is that already the first- 
order correction gives the exact form of the extremal fractional D-string solution (|35| ) . 
(j36j). Therefore, it is natural to expect that a similar expansion near the non-extremal 
D-string solution will capture the basic features of the non-extremal generalization of the 
extremal fractional D-string. 

More precisely, our starting point will be the well-known non-extremal ordinary D- 
string solution (f44|) with Q replaced by the effective D-string charge K*, so that we 
automatically match onto the near horizon asymptotics Perturbing in P 2 around the 
non-extremal D-string solution of charge A*, we will see that the 0(P 2 ) modification is 
already enough to match onto the extremal fractional D-string long-distance asymptotics. 
The small parameter governing this expansion is actually the dimensionless ratio A = 
P 2 K~ 1 a~ 1 ^ 3 , f\ i.e. for this method to work the horizon value of the effective D-string 
charge A* has to be sufficiently large. In view of the discussion in Section 5, this means 
that this perturbation theory is applicable for T 3> T c . Unlike ||, A here depends also 
on the non-extremality parameter a. 

It is useful to rescale the fields by appropriate powers of P 2 , setting 

K(u) = AT* + P 2 F(u) , $ n (u) = -2au/3 + P 2 <p(u) , w^u) = P 2 tOi(u) , (51) 

and 

y -> y + P 2 i , zn - z n + P 2 r] , (52) 

where y and z represent the pure D-string solution ([44]): e~ 6y = 3a" 1 sinh 2cm, e~ 4z ™ = 
^sinh2aw, and £ and r\ are corrections to it. To match onto the small u extremal 

3 Note that P ~ g s M and K <~ g s N where M and N are the numbers of fractional and regular 
D-strings respectively. 
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fractional D-string asymptotics (pE9"D, we require that 

c*(0) = m = m = , F-.-^_, 17 -^-p. (53) 

Now the system (^0|)— (25) takes the following explicit form: 

7£" + 0"-5O4e 12 ^ + O(P 2 ) = , (54) 

7wi' - <p" - 112e 12 ^i + 0{P 2 ) = , (55) 
2u% + 0" + 16e 12 ^ 2 + 0(P 2 ) = , (56) 
60" + e 42 "- 4 ^+ 2a «/ 3 (F' 2 - e ^- 4a "/ 3 ) + 0(P 2 ) = , (57) 

( e iz n -4y+2au/3 F y _ K ^BZn + (p2) = g (gg) 
8t/' - 4K*e 8Zn F - 16K 2 e 82 "r/ - e ^n-%+2au/3( F /2 + e 8y-4a«/3) + (p2) = q _ (gg) 

The constraint (^7| ) becomes 

16%'£' - 32*V + 4a0' - e ^-%+2a«/3( F /2 _ e 8^4a«/3) 

+ 8K 2 e 8z "r] + 2K*e 8Zn F - 84 • 12e 12 ^ + 0(P 2 ) = . (60) 

6.1 Leading-order solution for K 

Using the fact that K*e 8Zn = ^K' 1 ^ (g7]), we get from fl58|) that 

F' = e" 42 »- ta /3+%(C + 4K~ 1 <) . (61) 



For large u (near the horizon), we must have F' — > in order to satisfy (|45|). This along 
with the explicit form of z n ([44]) fixes the integration constant to be C = 2a/ K*. Hence 

/ 2a \ 2/3 
F' = -e~ 4au - ) , (62) 

V3(i- e - 4a «)y ' V ; 

and thus 

F =i^(J) (1-(1-^°") 1/3 ) • (63) 

As required by K(u) — K* + P 2 F{u) (|51~1), this expression satisfies F(w — > oo) = F* = 0. 
In other words, the fractional D-string charge K{u) — > in the large u limit, as desired. 
Notice that F' < so that, as advertised in Section 5, K(u) decreases as u increases. 
Moreover, if P 2 K~ 1 a~ 1 / 3 <C 1, the perturbation caused by F(u) is small for all values of 
u. Even at small u, F is well behaved 

3 /2\ 2 / 3 6 1 / 3 , 

F(u)= i^w -V" 1/3+o( " 4/3 )- < 64) 
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In addition, from these small u asymptotics, we find that at small u, we have recovered 
the extremal fractional D-string solution fl35|)! 

Thus already at the leading order this perturbation theory produces a solution with 
the correct extremal fractional D-string asymptotics. This remarkable fact strengthens 
our confidence that an exact solution interpolating between the extremal fractional D- 
string at small u and the regular non-extremal D-string horizon at large u indeed exists. 
Our perturbed solution should be a good approximation to it provided that P 2 K~ 1 <C 
a 1//3 . This limit corresponds to high Hawking temperatures, as we now show by matching 
(H) with (HD for small u. We find that 



3 • 6 1 / 3 



r -2 
Lp 



+ 



2\ 2 / 3 3 
3J 



(65) 



P 2 ' \3j 4oV3 

On the other hand, the Hawking temperature is determined in terms of the non-extremality 
a and the charge near the horizon K* by the usual D-string formula ( |? ) 

3 /3\ 1/6 
7T U, 



T 



(66) 



Comparing with the extremal fractional D-string solution, we expect that the critical 
temperature roughly corresponds to the value of a for which the naked singularity of 
the extremal fractional D-string solution and the horizon of the regular non-extremal 
fractional D-string coincide. In other words, a 1 ' 3 ~ Lp. Thus, the perturbation theory 
will correspond to high Hawking temperatures so long as a 1 / 3 ^> Lp. Comparing with 
(|65|), we recover the naive inequality needed for our perturbation theory to be valid, 
namely P 2 /K* < a 1/3 . To summarize, in our perturbative regime, T S> T c . 



6.2 Solutions for other fields 

We begin by solving for the correction to the dilaton 0. Using 
for the dilaton correction (1571) becomes 



1 

k1 



2a\ 5 / 3 
Ty (1 



-iau 



— e 



-Aau 



J 2/3 • 



Integrating once, and keeping in mind the boundary condition 



2a 



3fC 



4aK* 



2a \ 5 / 3 



and (H¥), the equation 



as u 



1 



-4au\l/3 



1 



(67) 

oo, we find 
(68) 



Hence, is a decreasing function of u. In other words, the string coupling decreases as 
we approach the horizon. To solve for 0, we may integrate once more: 



1 

8k* 



-3v 1/3 + I log — V m + y/3 tan 

2 1 — v 1 ' 6 



.! 2V 1 ' 3 + 1 vrv/3' 
V3 6~ 



(69) 
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From we have fixed the boundary condition such that 
to make the formula neater, we have introduced 



at u = 0. Moreover, 



-4cm 



and 



1 



2 . 2/3 



(70) 



(71) 



To show that <fi is well behaved at both small u and large u, we can write down some 
limiting expressions 



u — ► oo 



u — > 



+ id-) + o[(i-«) 2 ] 



1 



,4/3 



+0(t; 2 



where 



8fc* 3 + 



?r\/3 3 Q N 
^ + 2 l0g3 , 



(72) 
(73) 

(74) 



The relations describing the corrections to the other fields (|54])-(|56D, (|59|) are all 
ordinary, second order, linear differential equations. There exists a powerful technique, 
called the Lagrange method of variation of parameters, which is useful for dealing with 
this type of ODE. In particular, given an ODE of the form 



d 2 y ( ,dy 

— + p{x)— + q{x)y = g[x) 
ciar ax 



(75) 



and two linearly independent solutions y\ and y 2 to the corresponding homogenous equa- 
tion g[x) = 0, we can construct a general solution corresponding to the case g(x) ^ 0: 



Y = -Y 1 J ^-dx + Y 2 J ^-dx + c l Y l + c 2 Y 2 . 



Yig 



(76) 



where W = ^i^p — ^x~Y 2 is the Wronskian. (Linear independence corresponds to the 
fact that W^O.) 

As a first step, we cast the differential equations for £, and r], (|54l)-(|56|) and (|59|) , 
into the form ([75]) , i.e. 



(4 + 



64fl 2 e -4au 

9 (1 - e - 4au ) 2 
32a 2 e~ Aau 



2a 2 



-Aau 



21fc» (1 - e" 4aM ) 2 / 3 

,2 



a 



-4aw 



9 (1 



-4cmA2 



i0 2 



3K (1 - e- 4a «) 2 / 3 ' 



(77) 
(78) 



16 



v 



32a 2 e- Aau 



1 



a 
2k 



32a 2 e" 4au 
(1 

2 



-iau 



-4au\2 



12e 



-4au 



2a 2 

21K (1 - e- 4ffl «) 2 / 3 

lle -4a« 



-4aiA2 



-4em 



)5/3 



+ 



1-e 



-4au 



)5/3 



(79) 
(80) 



To analyze these differential equations, we again introduce the new radial variable v 
1 - e~ Aau . Then, Q77|)-(gg) become 



0>1 



^2 



4/9 



1 — v v 2 (l — V ) 

^ 2 2/9 

f 2 (1 — f ) 

2 



77 - 



7/ 



1 -v 

i 

l-v 
2 



^2 



f 2 (l — f) 
1 



1 - u 



-77 



168fc„ f 2 / 3 (l - u) 

1 1 

~48fc* v 2 / 3 (l-v) 

1 1 

_ 168^ v*/ 3 (l-v) 

12 11 



v 2 (l 



v 2 (l-v)'' 32k* 
where the dots denote d/dv. Now, the homogenous equation 

/ A 



,5/3 



/- 



1 — v v 2 (l — v) 



f = o 



(82) 
(83) 
(84) 

(85) 



is solved for generic A by /(f) = v v 2-^1 (^, v\ 2v\ v ), where 2 -Pi is the hypergeometric 
function and v{y — 1) = A. As it happens, A = 2 is a degenerate case where the 
solutions to the homogenous equation are elementary functions of v (namely, - — ^ an d 
—2 + log(l — v)). Perhaps surprisingly, we have been able to find analytic expressions 
for £ and 77 using these homogenous solutions. Let us start with £. In particular 



56k*v 



3v 1/3 (v - 4) + (u - 2) log 
+ 61 













+ b 2 { 









1 - 

-2 + 



'- + 2^(2 - v) tan" 1 



v/3 



+ 



log(l -v] 



(86) 



The boundary conditions for £ are that £ = at v = ( p3|) and d^/du — > as it — > 00 
This second boundary condition, in terms of the t> variable, means that d^/dv must 
be finite or zero in the limit v — > 1 because dv/du = 4a(l — u). As a result, we find that 



61 



7T 



>/3 



84*;* 



28*;* 



If) 



Next, we solve for 77: 
3 



64fc*t> 



67j 1/3 (f-2)-8 + (t»-2)log- 



f V3)3 2 7//3 

+ 2V3(2 - v) tan" 1 



v/3 
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+c 'U-2H(- 2+ — iog(i -^ 

Recall at the end of Section 5, we expected that the correction to the original dilaton 
$ should decrease as u — > oo, reflecting the corresponding decrease in K. Indeed, from 
fl6"9|) and (|88|), it is clear that this correction 30/2 — 2rj (see (0)) is indeed a decreasing 
function of it. The boundary conditions for rj are again that drj/du — > as u — > oo (45). 
However the small u boundary condition is merely the fact that 77 does not diverge at 
small u. The integration constants which satisfy these conditions are 

Cl = ^-(12-nV3); c 2 = . (89) 

To show that £ and 77 are well behaved at both large u and small u, we look at the 
asymptotics. First, we examine the long distance asymptotics (v — > or equivalently 
u — > 0) of £ and 77: 

^^ 4/3 -16k^ + ^ 7/3 )' ^ 

= - JL + + 9 ^4/3 _ 1 2 + 0( 7/3) (Q1) 

' 16/c* 6AK 896K G4h 

The first and second terms in the 77 expansion agree with the extremal fractional D-string 
solution, as one can see by comparing with ( pH ) and (EDI). 

Next we consider the expansions of £ and 77 near the horizon, u — > 00 or equivalently 
tj = 1: 

£ = £* + - -J—) (1 - t;) + 0[(1 - t;) 2 ] , (92) 



168A;, 
1 

32^ 



77 = 77* + ( 2r7» + ) (1 - v) + 0[(1 - vf] . ( 93 ) 



The horizon values 77* and are 



^ = T^r(- 15 + 7rv/ 3 + 91 °g 3 ) ' ( 94 ) 

168/c* 

77* = ^(-18 + 7rV3 + 91og3) . (95) 

We were careful to calculate £*, 77* and 0* because they show up as corrections to 
observables such as the Hawking temperature fl47|) and the horizon area. On the other 
hand, the volume of the compact space, this squashed Q 1 ' 1 ' 1 , does not depend on the 
squashing factors u\ and 0J2] the squashing factors cancel from the observables we are 
interested in calculating. As a result, and because the homogenous solutions to (|8l|) and 
(|S2j) are more complicated, we shall be more cavalier in our treatment of the asymptotics 
for u)\ and io-i- At small u, (|77|) and (|78|) become 

,,4/341/3 

Ux = PlU ^ + ^J_ u V3 logu + (u 7 ^) , (96) 
70 fc* 
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u 2 = p 2 u 1/3 + p 3 u 2/3 



,4/332/3 



-u 



4/3 



+ 0{u 



7/3 N 



(97) 



The Pi, i = 1,2,3, are undetermined constants of integration which can be determined 
numerically, keeping in mind the boundary conditions Wi(0) = and w'^oo) = 0. The 
large u asymptotics are 



1 



0>2 — ^2* 



168A;* 
1 



-iau 



+ 0(e 



-4au 



0(e 



-8au\ 



-8au\ 



(98) 



(99) 



We can now go a bit further than was done in the previous sections and determine the 
order P 2 corrections to the temperature and the entropy. The metric (HD can be recast 
into the form 



d* 2 



4a \ 3 / 4 
K~J 



e 3P 2 (v+H4) v -3/4 



v)dX 2 + dX\ 



+ 



4a\-V4 / 2a \V3 



KJ 



3 J 



1 

36 



fit \ u-i/12 [ J_ e P 2 m-v-<t>/4) 



dv 2 



v 2 {l — v) 
77 (p3|), and 



(dM 7 



(100) 



Using the large u asymptotics for £ fl92|), 77 Q93|) , and (f> (0), we obtain an explicit 
expression for the entropy per unit volume divided by the temperature squared: 



S 



VT 2 



a 



_^3/2 e 3P 2 (7^-2r ? ,-</»,/2) 



a 



K 3/2 ( ! , , D i n:i, 



8h 



(101) 



where a is a constant of order unity. In the last equality, we have used the values for 
(fH), 77* (|95|), and 0* (0). Note that the transcendental numbers log 3 and rr drop out 
of the linear combination of £*, 0*, and 77*. 
Using (|66D and flSlf) , one finds that 



p2 



1 - 



2 5 / 2 L 3 

3 l/2p T 



+ 



7T 



Hence, from (|101|) , 



S 



VT 2 



p3 
73" 



2 5 / 2 L 3 P 
7t3V 2 PT 



+ ... 



(102) 



(103) 



As expected, we find that the entropy ratio S/VT 2 increases toward a limiting value as 
T increases. The important point is that both the number of D-strings at the horizon 
and the entropy ratio S/VT 2 depend in the same way on the temperature. This picture 
is consistent with gauge theory where the number of D-strings should correspond roughly 
to the number of degrees of freedom in the theory. As the number of D-strings grows, so 
should the entropy. 
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7 Remarks 



We have presented, within the framework of perturbation theory in the number of frac- 
tional D-strings, a well behaved non-extremal fractional D-string solution. This finite 
temperature solution breaks supersymmetry but preserves the U(l) symmetry of the 
transverse conical space. It would be good if an additional KS type solution could be 
found for these fractional D-strings in a Q 1,1,1 background. In analogy with the KS so- 
lution for fractional D3-branes, we would expect the corresponding solution to involve 
blowing up the tip of the cone over Q 1 ' 1,1 in a way that keeps the five-cycle finite but 
allows the two cycle to become vanishingly small. In addition, we expect this solution 
would preserve supersymmetry but would break the U(l) symmetry of the Q 1 ' ' . If such 
a KS type solution is found, then we can use the finite temperature solution found here 
as evidence for a chiral symmetry breaking phase transition. We leave construction of 
such a deformation for future work. 

Another interesting direction to pursue is construction of a finite temperature solu- 
tion for the fractional D2-branes. In ||13|| , two examples of a KS type deformation of 
the fractional D2-branes were found. Each example involves a different six dimensional 
Einstein manifold - an S 2 bundle over either S 4 or CP 4 . It should be straightforward to 
introduce appropriate squashing functions to the cones over these two Einstein spaces, 
thereby producing an ansatz which admits well behaved fractional D2-branes at finite 
temperature, at least in perturbation theory. 

Finally, it would be extremely interesting to find a way to use the systems of differ- 



ential equations found here (|20l)-(p5|), ( [27]) and in |TT| to see what happens close to the 
expected phase transition. Unfortunately, to date we know of no such analytic solutions. 
Moreover, there are many integration constants involved, some of which are set in the 
IR, some in the UV, which makes any kind of numerical shooting algorithm a tedious 
prospect. Still, the authors of the current paper were mildly surprised that the set of 
differential equations @)~(25), (27) proved tractable in the first order in perturbation 
theory, and it seems likely that there are some more surprises waiting. 



Acknowledgments 

We are grateful to I. R. Klebanov and M. Rangamani for many useful discussions. The 
work of C. P. H. was supported in part by the DoD. P. O. is supported in part by an 
NSF Graduate Research Fellowship. 



References 



20 



J. Maldacena, "The Large N limit of superconformal field theories and supergravity," 
Adv. Theor. Math. Phys. 2 (1998) 231, frep-th/971120C| . 

S.S. Gubser, I.R. Klebanov, and A.M. Polyakov, "Gauge theory correlators from 



noncritical string theory," Phys. Lett. B428 (1998) 105, jhep-th/9802109 



E. Witten, "Anti-de Sitter space and holography," Adv. Theor. Math. Phys. 2 (1998) 
253, |hep-th/9802T5C. 



I. R. Klebanov and M. J. Strassler, "Supergravity and a Confining Gauge Theory: 
Duality Cascades and xSB-Resolution of Naked Singularities," JHEP 0008 (2000) 
052, frep-th/0007191] . 

S. S. Gubser, C. P. Herzog, I. R. Klebanov, and A. A. Tseytlin, "Restoration of 
Chiral Symmetry: A Supergravity Perspective," [hep-tri/0102172[ 



I. R. Klebanov and N. Nekrasov, "Gravity Duals of Fractional Branes and Logarith- 



mic RG Flow," Nucl. Phys. B574 (2000) 263, [hep-th/99lT096 



I.R. Klebanov and E. Witten, "Superconformal field theory on threebranes at a 
Calabi-Yau singularity," Nucl. Phys. B536 (1998) 199, [hep-th/ 9807080. 



D. Morrison and R. Plesser, "Non-Spherical Horizons, I," Adv. Theor. Math. Phys. 



3 (1999) 1, |hep-th/98 10201 



I. R. Klebanov and A. A. Tseytlin, "Gravity Duals of Supersymmetric SU(iV) x 
SU(iV + M) Gauge Theories," Nucl. Phys. B578 (2000) 123, |hep-th/0002T5S| . 

A. Buchel, "Finite temperature resolution of the Klebanov- Tseytlin singularity," 
piep-th/0011146 . 



A. Buchel, C. P. Herzog, I. R. Klebanov, L. Pando Zayas, and A. A. Tseytlin, "Non- 



Extremal Gravity Duals for Fractional D3-Branes on the Conifold," [hep-th/0102105 



C. P. Herzog and I. R. Klebanov, "Supergravity Duals of Fractional Branes in Various 
Dimensions," |iep-th/0101020 . 



M. Cvetic, G. W. Gibbons, H. Lii, and C. N. Pope, "Supersymmetric Non-Singular 
Fractional D2-branes and NS-NS 2-branes," [hep-th/0101096| . 

L. Castellani, L. J. Romans, and N. P. Warner, "A Classification of Compactifying 
Solutions for d=ll Supergravity," Nucl. Phys. B241 (1984) 429. 

A. Buchel and A. Frey, "Comments on Supergravity Dual of Pure N=l Super Yang- 
Mills Theory with Unbroken Supersymmetry," [hep-th/P 103022 . 



21 



[16] J. Maldacena and C. Nunez, "Towards the Large N Limit of Pure N=l Super Yang- 
Mills," Phys. Rev. Lett. 86 (2001) 588, frep-th/0008001 



[17] D. J. Gross and E. Witten, "Possible Third Order Phase Transition in the Large N 
Lattice Gauge Theory," Phys. Rev. D. 21 (1980) 446. 

[18] M. Cvetic, G. W. Gibbons, H. Lii, and C. N. Pope, "Ricci-flat Metrics, Harmonic 
Forms and Brane Resolutions," [hep-th/0012011| . 

[19] G. T. Horowitz and A. Strominger, "Black strings and P-branes," Nucl. Phys. B360 
(1991) 197. 

[20] N. Itzhaki, J. M. Maldacena, J. Sonnenschein, and S. Yankielowicz, "Supergravity 
and the Large N Limit of Theories with Sixteen Supercharge," Phys. Rev. D 58 
(1998) 046004, [hep-th/9802042 . 



22 



